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Summary. — We search for new charmless decays of neutral fe-hadrons to pairs 
of charged hadrons with the upgraded Collider Detector at the Fermilab Tevatron. 
Using a data sample corresponding to lfb - of integrated luminosity, we report the 
first observation of the B® — > K~tt + decay, with a significance of 8.2a, and measure 
B(B° -> K-n + ) = (5.0 ± 0.7 (stat.) ± 0.8 (syst.)) x 10 -6 . We also report the first 
observation of charmless fe-baryon decays in the channels Aj, — > piv~ and A° — > pK~ 
with significances of 6.0<r and 11.5<r respectively, and we measure £>(A° — > ptv~) = 
(3.5 ± 0.6 (stat.) ± 0.9 (syst.)) x 10" 6 and B(A° b pK~) = (5.6 ± 0.8 (stat.) ± 
1.5 (syst.)) x 10" 6 . No evidence is found for the decays B° — > K + K~ and B° — > 
tv + tt~ , and we set an improved upper limit B(B® — } 7r + 7r _ ) < 1.2 x 10~ 6 at the 90% 



confidence level. All quoted branching fractions are measured using B(B° 
as a reference. 
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1. Introduction 

Non-lcptonic two-body charmless decays of neutral b hadrons (B° -> h+h S° -> 
h + h~ and A° — > p/i - , where ft, is a charged pion or kaon) are very interesting for the 
understanding of flavor physics and CP violation mechanism in the 6-hadron sector. Their 
rich phenomenology offers several opportunities to explore and constrain the parameters 
of the quark-mixing matrix (i.e. Cabibbo-Kobayashi-Maskawa, CKM). These processes 
allow to access the phase of the V u b element of the CKM matrix (7 angle) , and to test the 
reliability of the Standard Model (SM) and hadronic calculations. The presence of New 
Physics can be revealed by its impact on their decay amplitudes, where new particles 
may enter in penguin diagrams. The B° — > K + tt~ is the first process involving the b 
quark where direct CP violation has been observed. 

The measurements obtained at e + e~ colliders (ARGUS, CLEO, LEP, and more re- 
cently, BaBar and Belle experiments) already provided a wealth of results for B° and B + 
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mesons. The upgraded Collider Detector at the Fermilab Tevatron (CDF II), with its 
large production of 6-hadrons is in principle an ideal environment for studying these rare 
modes. In addition to providing further large samples of B° and B + mesons in a different 
experimental environment, it provides the exciting opportunity of studying the charmless 
decays of other 6-hadrons that are unaccessible (or much less accessible) in other exper- 
iments. A variety of techniques have been proposed to constrain the CKM parameters 
or probe effects of New Physics [1, 2, 3, 4] exploiting a combination of observables from 
B° s and B° , B+ mesons. 

The B° s — > K~tt + decay offers several and interesting strategies to extract useful 
information from the comparison between its observables and those of its U-spin related 
partner B a — > K + n~ [3, 4]. By combining the information of rates and direct CP 
asymmetries of U-spin- -related decays B° -> K+tt~ and B° s -> K~ir + [3, 4] allows a 
stringent test of the Standard Model origin of the 0(10%) direct CP asymmetry observed 
in B° — > K + tt~ [8], which is not matched by a similar effect in the B + — > K + n° decay, 
which differs only by the spectator quark. This raised discussions about a possible 
exotic source for the CP violation in the B° — > K + n~ decay [9, 10, 11]. Any significant 
disagreement between the measured partial rate asymmetries of strange and non-strange 
6-meson Kit decays should be strong indication of New Physics. 

The B° s — »■ K~ir + is still unobserved and the current experimental upper limit 
B(B° S -> K-TT+) < 5.6 x 10~ 6 @ 90% CL [5] from CDF is very close to (sometimes 
lower than) the current theoretical expectations [12, 13, 14, 15]. The comparison of this 
branching fraction, sensitive to CKM angle values of a and 7 [14], with theoretical predic- 
tions provides valuable information for tuning the phenomenological models of hadronic 
B9 s decays and for optimizing the choice of their input parameters. Therefore, in this 
context the measurement of the decay rate of the B® — > K~n + and the measurement of 
its direct CP asymmetry becomes crucial. 

The amplitudes of penguin-annihilation and exchange diagrams, in which all initial- 
state quarks undergo a transition, are difficult to predict with current phenomenological 
models. In general they may carry different CP-violating and CP-conserving phases with 
respect to the leading processes, thereby influencing the determination of CKM-related 
parameters. The B° — > K + K~ and B° — »■ ir + -K~ decays proceed only through these kinds 
of diagrams. A simultaneous measurement of their decay rates (or improved constraints 
on them) would provide valuable estimates of the magnitude of these contributions [16]. 

Simultaneous measurements of B^ — > h + h ~ observables, in most cases, exploit the 
U-spin symmetries to partially cancel out or constrain hadronic uncertainties and probe 
the electroweak and QCD structure. U-spin symmetry is not exactly conserved in the 
Standard Model and the magnitude of its violation is not precisely known but most 
authors estimate a 0(10%) effect. The B?^ — > h + h ~ system is a privileged laboratory 
since it offers the simultaneous opportunities of using U-spin assumptions and, at the 
same time, of checking their validity by measuring the symmetry breaking-size, from the 
interplaying of several U-spin-related observables. 

Two-body charmless decays are also expected from bottom baryons. The modes 
A° — > pK ~ and A° — > pir~ are predicted to have measurable branching fractions, of 
order 10~ 6 [17], and, in addition to the interest in their observation, must be considered 
as a possible background to the rare B° and B° modes being investigated. 

In this Letter we report the results of a search for rare decays of neutral bottom 
hadrons into a pair of charged charmless hadrons (p, K or tt), performed in lfb -1 of pp 
collisions at y/s = 1.96 TeV, collected by the upgraded Collider Detector (CDF II) at the 
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Fcrmilab Tevatron. We report the first observation of modes B° — > K~tt + , A" — > pK~ , 
and A° — > p7r~, and measure their relative branching fractions. This is a short overview 
of the work documented in Ref. [18] and the results were published in [19]. 

Throughout this paper, C-conjugate modes are implied and branching fractions indi- 
cate CP-averages unless otherwise stated. 

2. CDF II detector 

The CDF II detector [20, 21], in operation since 2001, is an azimuthally and forward- 
backward symmetric apparatus designed to study pp collisions at the Tevatron. It is a 
general purpose solenoidal detector which combines precision charged particle tracking 
with fast projective calorimetry and fine grained muon detection. Tracking systems 
are contained in a superconducting solenoid, 1.5 m in radius and 4.8 m in length, which 
generates a 1.4 T magnetic field parallel to the beam axis. Calorimetry and muon systems 
are all outside the solenoid. The main features of the detector systems are summarized 
below. 

The tracking system consists of a silicon microstrip system [22] and of an open-cell 
wire drift chamber [23] that surrounds the silicon. The silicon microstrip detector con- 
sists of seven layers (eight layers for 1.0 < |ry| < 2.0) in a barrel geometry that extends 
from a radius of r = 1.5 cm from the beam line to r — 28 cm. The layer closest 
to the beam pipe is a radiation-hard, single sided detector called Layer00 which em- 
ploys sensors supporting high-bias voltages. This enables signal-to- noise performance 
even after extreme radiation doses. The remaining seven layers are radiation-hard, dou- 
ble sided detectors. The first five layers after Layer00 comprise the Silicon VerteX 
II (SVXII) detector and the two outer layers comprise the Intermediate Silicon Layer 
(ISL) system. This entire system allows track reconstruction in three dimensions. The 
impact parameter resolution of the combination of SVXII and ISL is about 48 ^m in- 
cluding a 30 ^m contribution from the beamline for tracks with transverse momentum 
of 2 GcV/c. The z resolution of the SVXII and ISL is 70 (j,m. The 3.1 m long cylin- 
drical drift chamber (COT) covers the radial range from 40 to 137 cm and provides 96 
measurement layers, organized into alternating axial and ±2° stereo superlayers. The 
COT provides coverage for |r;| <1. The hit position resolution is approximately 140 /zm 
and the momentum resolution <j Pt /pt — -15% pt / (GeV / c) . This corresponds to an 
observed mass-widths of about 14 MeV/c 2 for the J/tjj — » n + ^~ decays, and of about 8 
MeV/c 2 for the D° — > K~ir + decays. The specific energy loss by ionization (dE/dx) of 
charged particles in the COT can be measured from the amount of charge collected by 
each wire. This yields a nearly-constant separation of 1.5 standard deviations between 
pions and kaons over the range 2 < px < 10 GeV/c. 

A Time-of-Flight (TOF) detector [24], based on plastic scintillators and fine-mesh 
photomultiplicrs is installed in a few centimeters clearance just outside the COT. The 
TOF resolution is « 100 ps and it provides at least two standard deviation separation 
between K ± and ^ for momenta p < 1.6 GcV/c. 

Segmented electromagnetic and hadronic sampling calorimeters surround the tracking 
system and measure the energy flow of interacting particles in the pseudo-rapidity range 
\r]\ < 3.64. The central calorimeters (and the endwall hadronic calorimeter) cover the 
pseudorapidity range \t]\ < 1.1(1.3). The central electromagnetic calorimeter [25] (CEM) 
uses lead sheets interspersed with polystyrene scintillator as the active medium and em- 
ploys phototube readout. Its energy resolution is 13.5%/v / ^t©2%. The central hadronic 
calorimeter [26] (CHA) uses steel absorber interspersed with acrylic scintillator as the 
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active medium. Its energy resolution is 75% /\/~Et ® 3%. The plug calorimeters cover 
the pseudorapidity region 1.1 < \r/\ < 3.64. They are sampling scintillator calorimeters 
which are read out with plastic fibers and phototubes. The energy resolution of the plug 
electromagnetic calorimeter [27] is 16%/y/E © 1%. The energy resolution of the plug 
hadronic calorimeter is 74%/ VE ® 4%. 

The muon system resides beyond the calorimetry. Four layers of planar drift chambers 
(CMU) detect muons with p T > 1.4 GeV/c which penetrate the five absorption lengths of 
calorimeter steel. An additional four layers of planar drift chambers (CMP) instrument 
0.6 m of steel outside the magnet return yoke and detect muons with px > 2.0 GeV/c. 
The CMU and CMP chambers each provide coverage in the pseudo-rapidity range < 
0.6. The Intermediate MUon detectors (IMU) are covering the region 1.0 < \t]\ <1.5. 

The beam luminosity is determined by using gas Cherenkov counters located in the 
3.7 < \r]\ < 4.7 region which measure the average number of inelastic pp collisions per 
bunch crossing [28]. 

The trigger and data acquisition systems are designed to accommodate the high 
rates and large data volume of Run II. Based on preliminary information from tracking, 
calorimetry, and muon systems, the output of the first level of the trigger is used to limit 
the rate for accepted events to s=s 18 kHz at the luminosity range of 3-7 10 31 cm _2 s _1 . 
At the next trigger stage, with more refined information and additional tracking infor- 
mation from the silicon detector, the rate is reduced further to w 300 Hz. The third and 
final level of the trigger, with access to the complete event information, uses software 
algorithms and a computing farm, and reduces the output rate to w 75 Hz, which is 
written to permanent storage. 

The only physics objects used in this analysis are the tracks, then just tracking system 
has been used. 

3. Data sample 

We analysed an integrated luminosity J Cdt ~ 1 fb _1 sample of pairs of oppositely- 
charged particles with pr > 2 GeV/c and pn +PT2 > 5.5 GeV/c, used to form 6-hadron 
candidates. The trigger required also a transverse opening angle 20° < A<f) < 135° 
between the two tracks, to reject background from particle pairs within the same jet and 
from back-to-back jets. In addition, both charged particles were required to originate 
from a displaced vertex with a large impact parameter d (100 /im < d < 1 mm), while 
the 6-hadron candidate was required to be produced in the primary pp interaction (ds < 
140 /im) and to have travelled a transverse distance Lt > 200 ttm. 

The offline selection is based on a more accurate determination of the same quantities 
used in the trigger, with the addition of two further observables: the isolation (Ib) of 
the candidate [29], and the quality of the three-dimensional fit (x 2 with 1 d.o.f.) of the 
decay vertex of candidate. Requiring a large value of Ib reduces the background from 
light-quark jets, and a low \ 2 reduces the background from decays of different long-lived 
particles within the event, owing to the good resolution of the SVX detector in the z 
direction. 

In the offline analysis, an unbiased optimization procedure determined a tightened 
selection on track-pairs fit to a common decay vertex. We chose the selection cuts mini- 
mizing directly the expected uncertainty of the physics observables to be measured. The 
selection is optimized for detection of the B® — > K ~it + mode. Maximal sensitivity for 
both discovery and limit setting is achieved with a single choice of selection require- 
ments [30] by minimizing the variance of the estimate of the branching fraction in the 
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absence of signal [18]. The variance is evaluated by performing the full measurement 
procedure on simulated samples containing background and all signals from the known 
modes, but no B® — > K~ir + signal. The variance has been parameterized with analytical 
functions of the signal yield (S) and background level (B), and the free parameters deter- 
mined from analysis of pseudo-experiments reproducing the experimental circumstance 
of data. For each set of cuts, S was estimated from Monte Carlo simulation and nor- 
malized to the yield observed in data after the trigger selection, and B was extrapolated 
from the sidebands of the 7T7r-mass distribution in data. This procedure yields the final 
selection: Ib > 0.525, x 2 < 5, d > 120 ^m, d,B < 60 ^m, and Lt > 350 /xm. 

The resulting invariant-7T7r-mass distribution (see Fig. 1) shows a clean signal of 
B9-. — ¥ h + h ~ decays. In spite of a good mass resolution 22MeV/c 2 ), the various 
signal decay modes overlap into an unresolved mass peak. 




Invariant 7nt-mass [GeV/c 2 ] Invariant 7nt-mass [GeV/c 2 ] 

Fig. 1. - Invariant 7T7r-mass distribution of reconstructed candidates. The charged pion mass 
is assigned to both tracks. The total projection and projections of each signal and background 
component of the likelihood fit are overlaid on the data distribution. Signals and multi-body B 
background components are shown stacked on the combinatorial background component. Linear 
scale (a), logarithmic scale (b). 



4. Fit of composition 

The resolution in invariant mass and in particle identification is not sufficient for sep- 
arating the individual decay modes on an event- by- event basis, therefore we performed 
an unbinned maximum likelihood fit, combining kinematic and particle identification in- 
formation to statistically determine both the contribution of each mode, and the relative 
contributions to the CP asymmetries. For the kinematic portion, we used three loosely 
correlated observables to summarize the information carried by all possible values of in- 
variant mass of the 6-hadron candidate, resulting from different mass assignments to the 
two outgoing particles [31]. They are: (a) the invariant 7T7r-mass m^n calculated with the 
charged pion mass assignment to both particles; (b) the signed momentum imbalance 
a = (1 — p\/p2)qi 1 where p\ (j>2) is the lower (higher) of the particle momenta, and qi is 
the sign of the charge of the particle of momentum pi ; (c) the scalar sum of the particle 
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momenta ptot = Pi + Pi- Using these three variables, the mass of any particular mode 
m mim2 (m K7r ,m nK , m KK ,m pv , m^p, m pK ,m Kp ) can be written as: 



(1) 



2mJ + (m 1 + m 2 ) 



(2) 



Pi = 



2 - a 



■Ptot , Pi = 



2- la 



■Ptot, 



where m\ (77J2) is the mass of the lower (higher) momentum particle. For simplicity, 
Eq. (1) is written as a function of p\ and p 2 , but in the likelihood it was used as a 
function of a and p tot . The simulated average values of function of a for the 

twelve — > h + h~ and A" — > p/i~ modes are shown in Fig. 2 and Fig. 3. Particle 
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Fig. 2. - Average m^vr versus a for simulated samples of B°, and A° decay modes. 



identification (PID) information is summarized by a single observable kaonness ni(i) for 
track 1(2), defined as 

dE/dx^) - dE/dx 1( 2){n) 
Kl(2) " dE/dx 1(2) {K)-dE/dx m (ir)' 

where dE/dxi^) {^) and dE/dxi^)(E) are the expected dE/dxip) depositions for those 
particle assignments. With the chosen observables, the likelihood contribution of the i th 
event is written as: 

(3) d = (l- f b ) /i4 in 4 ID + h (/A£ A in ^A D + (1 - /A)£ C in £c D ) 

j 

where: 

(4) C) m = Rj(m nv \a,ptot)Pj(a,Ptot), 



(5) 



£ A m = A(m 7r7r ;c 2 ,m )PA(a,Ptot), 



OBSERVATION OF NEW CHARMLESS DECAYS OF BOTTOM HADRONS 



rr- 5.6: 

? 5.55 r B° 



R) 5-45 

E E „ 



! 5.3 
i 

I 5.25 
5.2 
5.15 



5.1 



■0.6 -0 4 -0.2 -0 0.2 0.4 0.6 0.8 1 

a = (1 - p/p 2 ) x q, 



p-i 5.6: 

S 5.55^ B° 

0) - 

a 5.5^ 

tf> 

I « 

N 

_ 5.35 
C 

.2 5.3 

a 

£ 5.25 
5.2 
5.15 



5.1 



-0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 1 

a = (1 - P/P 2 ) x q i 



F-i 5.6 p 

3? : =o 

> 5.55 ,B 
£ 5.5 r 

(0 

SS"5r 

I M 

h; 

„ 5.35 
C 

5 5.3 
CO 

> 5.25 

5.2 
5.15 



5.1 



1 -0.8 -0.6 -0.4 -0.2 -0 0.2 4 0.6 0.8 1 

a = (1 - p/p 2 ) x q t 



T 56 : 
> 5.55: 



B° -» K + K" 



3 5-45 r 

| 5 - 4 r 

5.35 r 
3 5.3^ 

3 

i 5.25 r 
5.2 r 
5.15 r 




0.2 0.4 0.6 0.8 1 

a = (1 - p,/p 2 ) x q ( 



> 5.55 B* 

a 5.5 
tf) 

IS MS 

I " i 

_ 5.35 

C 

.2 5.3 
CO 

> 5.25 

5.2 



-0.6 -0.4 -0 2 -0 0.2 0.4 0.6 0.8 1 

a = (1 - p/p 2 ) x q| 



?5.55rB°^KV 



«5.45^ 

E C A- 



-1 -0.8 -0.6 -0.4 -0.2 -0 



0.2 0.4 
X = (1- 



0.6 0.8 1 

P/P 2 )xq, 



rr^ 5.6 — 

> 5.55 - I 
tf) 

fl 5.5^ 
tf) 

S M5r 
1 M 

^ 5.35 r 
£ 

.2 5.3 r 
to 

> 5.25 r 

5.2 r 
5.15 r 



5.1? 



-0.6 -0.4 -0.2 -0 



0.2 0.4 0.6 0.8 1 

x = (1 - p/p 2 ) x q. 



> 5.55 r B" 
eo s 

s 5.5 

I MS 

1 54 

_ 5.35 
C 

.2 5.3 
CO 

> 5.25 

5.2 r 
5.15 r 
5.1^ 



>K + K" 



-0.6 -0.4 -0.2 -0 



0.2 0.4 0.6 0.8 1 

i = 0- p/p 2 ) x q, 



S5.45, 

| 5.4, 

^ 5.35 r 
C 

(0 

£ 5.25 r 
5.2 r 
5.15 r 





AS >P 









5.1 



1 -0.8 -0.6 -0.4 -0.2 -0 



0.2 0.4 
X = (1- 



0.6 0.8 1 

P,'P 2 )xq, 



> 5.S5 A b 
2. 5.5 

tf) 

<S M5 
I " 

R 

_ 5.35 
C 

.2 5.3 
CO 

g 5.25 
5.2 









■ 







-0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 1 

a = (1 - P,/p 2 ) x q ( 



rr 5.6: 

> 5.55 : 
0) 

a 5.5^ 

in 

8 5.45^ 



A° — * pK' 



■0.6 0.4 0.2 -0 0.2 0.4 0.6 0.8 1 
« = (1 • p/p ; ) x q| 



»5.45^ 



,pK* 



■1 0.8 0.6 0.4 0.2 



0.2 0.4 0.6 0.8 1 

» = (1 - p,'p 2 ) x q, 



Fig. 3. - m^-Tr versus a for simulated samples of B? s -, h + h and A° — > ph decay modes. 
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(8) C^ = ]T w^ C) w^F lm ( Ku K 2 \a,p tot ). 

l,m—e,7T,K,p 

The various terms of the likelihood functions are described below. 

The index j runs over the twelve distinguishable £?9n — > h + h~ and A" — > modes, 
and fj are their fractions to be determined by the fit, together with the total back- 
ground fraction The background is composed of two different kinds: combinatorial 
background and partially-reconstructed heavy flavor decays. The combinatorial back- 
ground is composed of random pairs of charged particle, displaced from the beam-line, 
accidentally satisfying the selection requirements, while the latter, referred as "physics" 
background, is composed of multi-body &-hadron decays (i.e. £?9n — > pir/pK) in which 
only two tracks are reconstructed. The indices A(C) label the physics (combinatorial) 
background quantities. The fraction of the physics background is given by /a and it is 
a free parameter in the fit. 

Each likelihood term, both for signals and backgrounds, is factorized into three dif- 
ferent contributions: a) the conditional probability distribution of the invariant mass 
m^ given a and ptot (for the background m V7r is assumed to be independent of mo- 
mentum), b) the joint conditional probability of PID variables m, n 2 given a, p to t f° r 
a determined particles hypothesis, j in the case of signals (Fj) and l,m in the case of 
background (F^ m ), and c) the joint probability distribution of momentum variables a 
and p tot (P,-(a,c))- 

If TV (rrij ) is the mass resolution function of each mode j when the correct mass is 
assigned to both tracks, we can use Eq. (1) to change variable rrij — > m W7r and to write 
the density probability function for each j th decays mode as function of m„ given a and 
Ptot- In fact: 

din 

(9) Rj(mj) = Tlj(mj(m vw )) ■ ^ = Rj(m vn \a,p tot ). 

The functional form of the mass resolution function IZj (rrij ) was parameterized using the 
detailed detector simulation. To take into account non-Gaussian tails due to the emission 
of photons in the final state, we included soft photon emission in the simulation, using 
recent QED calculations [32]. The quality of the mass resolution model was verified 
by comparison data and simulation with about 1.5 x 10 6 tagged D° — > K~n + decays 
reconstructed using the chain D* + — > D°ir + — )• [K~tt + ]tt + (see Fig. 4 and Fig. 5). 
The mass line-shape of the D° — > K~ir + was fitted by fixing the signal shape from the 
model, and allowing to vary only the background function. Good agreement was obtained 
between data and simulation. In Eq. (4), the nominal B°, B® and A° masses measured 
by CDF [33] were used to reduce the systematic uncertainties related to the knowledge 
of the global mass scale. 

The mass distribution of the physics background is parameterized with an "Argus 
function", defined by the notation A(m 7r7r ; c 2 , mo) [34], convoluted with a Gaussian dis- 
tribution centered at zero with a width, in this case, equal to the mass resolution, while 
the combinatorial background with an exponential function. The background mass dis- 
tribution was determined in the fit by varying the parameters c\, c 2 and mo in Eq. (5,6). 
The function Pj(A.c) i a ,Ptot ) was parameterized by a product of polynomial and exponen- 
tial functions fitted to Monte Carlo samples produced by a detailed detector simulation 
for each mode j, instead for the background terms was obtained from the mass sidebands 
of data [18]. 
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Fig. 4. - Invariant-if7r mass distribution for tagged D —> K~tt + decays from D* + — ► D % + — > 
[K~ ir + ]-K + . A verification of the mass line shape is superimposed, by performing a 1-D binned 
fit where the signal mass line shape is completely fixed from the model (see text). Linear scale 
(a), logarithmic scale (b). 



The same data sample of 1.5 x 10 6 D* + -» D°ir + — s- [_ftT~7r + ]-7r + decays used to test 
mass resolution model, where the D° decay products are identified by the charge of the 
D* + pion, was used to calibrate the dE/dx response over the tracking volume and over 
time, and to determine the Fj^ m )(ni, K2\a,ptot) functions in Eq. (7,8). In a > 95% pure 
D° sample, we obtained approximately 1.5er separation between kaons and pions for par- 
ticles with momentum larger than 2 GeV/c (see Fig. 6), corresponding to an uncertainty 




Fig. 5. - Invariant-7f7r mass distribution of simulated B° —¥ K + ir~ decays. The mass template 
is superimposed. Linear scale (a), logarithmic scale (b). Similar templates for all B? e \ — >• h + h ~ 
and — > ph~ decay modes. 
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Fig. 6. - Distribution of dE/dx around the average pion response for negatively- (a) and 
positively- (b) charged particles. Pions (continuous line) and kaons (dashed line) from D° — > 
K~ir decays. 



on the measured fraction of each class of particles that is just 1.7 times worse than the 
uncertainty attainable with ideal separation between two classes of events completely 
disentangled. The effective separation among final states consisting in particle pairs, 
like in our case (between 7r + 7r~ and K + K~, between ir + K~ and K + tt~) corresponds 
to 1.5a ■ \/2 ~ 2.1c, as shown in Fig. 7. This achievement is particularly crucial in 
separating those signal decay modes in which the kinematics does not sufficiently help. 



3 if^Kn 



-ir 

-2 7 
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, ■ ■■■■■■■■ ■ ■ ■ . . 
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■ ■ ■ ■■■■■■ -J □ □ □ □ B 

JJJJ-. DDDDBI 



3 B°/B° K + K 
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K, 



Fig. 7. - Probability density function J Fj(n\, K,2\a, Ptot)Pj (a, ptot)dadptot for the some signal 
decays in the space a > 0. To obtain the distributions for a < it is sufficient to invert «i <H> k-2 
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Table I. - Yields of signals returned from fits. For rare unobserved modes significance is quoted. 
The first quoted uncertainty is statistical, the second is systematic. 



Mode 




Significance 


B° -> K--K+ 


230 ± 34 ± 16 


8.2a 


B° -> TV + TV- 


26 ± 16 ± 14 


< 3(7 


B° -> K+K~ 


61 ± 25 ± 35 


< 3(7 


A° b -> pK~ 


156 ± 20 ± 11 


11.5(7 




110 ± 18 ± 16 


6.0a 



For example the kinematic separation power between the B° — > 7r + 7r~ and — > K + K~ 
modes is almost null, as shown in Fig. 2 and Fig. 3, while the dE/dx power separation 
is maximum, about 2.1a, as shown in Fig. 7. 

The dE/dx response of protons was determined from a sample of 124,000 A — > p7r~ 
decays, where the kinematics and the momentum threshold of the trigger allow unam- 
biguous identification of the decay products [18]. 

The PID background term in Eq. (8) is similar to the signal terms, but allows for 
independent pion, kaon, proton, and electron components, which are free to vary inde- 
pendently for physics (combinatorial) background. In Eq. (8) the indices I and m run 
over the four possible particles e, ir, K, p and the fractions of different kind of particles 
wf^ E \wm^ are free parameters in the fit. Muons are indistinguishable from pions with 
the available dE/dx resolution. 

From the signal fractions returned by the likelihood fits we calculate the signal yields 
shown in Table I. The significance of rare unobserved signals is evaluated as the ratio 
of the yield observed in data, and its total uncertainty (statistical and systematic) as 
determined from a simulation where the size of that signal is set to zero. This evaluation 
assumes a Gaussian distribution of yield estimates, supported by the results obtained 
from repeated fits to simulated samples. This procedure yields a more accurate measure 
of significance with respect to the purely statistical estimate obtained from y/— 2Aln(£). 
Significant signals are seen for B a — > tt + it~ , B — > K + n~ , and B® — > K + K~ , previously 
observed by CDF [5]. We obtain significant signals for the B® — »■ K~n + mode (8.2<r), 
and for the A" — > pir~ (6.0<r) and A" — > pK~ (11.5er) modes. Figure 8 shows relative 
likelihood distributions for these modes. No evidence is found for the modes B® — > tt + it~ 
or B a — > K + K~, in agreement with expectations of significantly smaller branching 
fractions. 

To avoid large uncertainties associated with production cross sections and absolute 
reconstruction efficiency, we measure all branching fractions relative to the B° -4- K + ir~ 
mode. Frequentist upper limits [35] at the 90% C.L. are quoted for the unseen modes. 
For the measurement of A£ branching fractions, the additional requirement j»t(A°) > 
6 GeV/c was applied to allow easy comparison with other A^J measurements at the Teva- 
tron, which are only available above this threshold [36, 37]. This additional requirement 
lowers the A° yields by about 20%. 

5. Acceptance corrections and systematics 

To convert the yields returned from the fit into relative branching fractions measure- 
ments, we applied corrections for efficiencies of trigger and offline selection requirements 
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Fig. 8. - Distribution of the relative signal likelihood, Cs/{£-s + £ thcr), in the region 5.1 < 
m-Kir < 5.6 GeV/c 2 . For each event, Cs is the likelihood for the B° — > K~n + (a), A° — > 
pK~ (b), or — > p-K~ (c) signal hypotheses, and £ thcr is the likelihood for everything but 
the chosen signal, i.e. the weighted combination of all other components according to their 
measured fractions. Points with error bars show the distributions of data and histograms show 
the distributions predicted from the measured fractions. 



for different decay modes. The relative efficiency corrections between modes do not ex- 
ceed 8% for the measurements of 6-mesons and 40% for A° branching fractions. Most 
corrections were determined from the detailed detector simulation, with some exceptions 
which were measured using data. A momentum-averaged relative isolation efficiency 
between and B° of 1.00 ± 0.03 was determined from fully-reconstructed samples of 
B° s — > J/ij^cf) an d B°-> J/ipK*° [18]. The lower specific ionization of kaons with respect 
to pions in the drift chamber is responsible for a ~ 5% lower efficiency to reconstruct a 
kaon. This effect was measured in a sample of D + — > K~it + it + decays triggered on two 
tracks, using the unbiased third track [38]. 

The B° s — > n + ir~ modes required a special treatment, since it contains a superposition 
of the flavor eigenstates of the B Q S . Their time evolution might differ from the one of 
the flavor-specific modes if the width difference AT S between the B® mass eigenstates 
is significant. The current result was derived under the assumption that both modes 
are dominated by the short-lived B° s component, that Y s = Yd, and Ar s ./r s = 0.12 ± 
0.06 [39, 40]. The latter uncertainty is included in estimating the overall systematic 
uncertainty. 

The dominant contributions to the systematic uncertainty are the uncertainty on 
the combinatorial background model and the uncertainty on the dE/dx calibration and 
parameterization. Smaller systematic uncertainties are assigned for trigger efficiencies, 
physics background shape, kinematics, B meson masses and lifetimes. 

6. Results 

The final results on branching fractions are listed in Table II, where fd, f s and /a 
indicate the production fractions respectively of B°, B® and A° from fragmentation of a 
b quark in pp collisions. An upper limit is also quoted for modes in which no significant 
signal is observed [35] . We also list absolute results obtained by normalizing the data to 
the world-average of B(B° — > K + ir~) [43, 44]. The contributions from the likelihood fit 
for each decay mode are shown in Fig. 1. 

The CP-averaged branching fraction of the newly observed mode B® — > K~ir + is 



OBSERVATION OF NEW CHARMLESS DECAYS OF BOTTOM HADRONS 



13 



Table II. - Measured relative branching fractions of rare modes. The ratio fh/fd is pr~ 
dependent [37], and is defined here as: fh/fd = cr(pp — > A°X;pt > 6 GeV/c, \r/\ < l)/a(pp — > 
B°X;pr > 6 GeV/c, \r/\ < 1). Absolute branching fractions were derived by normalizing to the 
current world-average value B{B° -> K+iv~) = (19.4±0.6) x 10~ 6 , and assuming the average val- 



ues at high energy for the production fractions: f s /fd = 0.276±0.034, ana 
[43]. The first quoted uncertainty is statistical, the second is systematic. 


',f\/U = 0.230 ±0.052 


Mode 




Relative B 


Absolute S(10 -6 ) 


B° -> K-K+ 
B° ->■ n+n- 

B° ->■ K+K~ 

A° -»• pK- 
A% ->■ P 7Y- 


/, B{B a s ^K-TT+) 


0.071 ± 0.010 ± 0.007 
0.007 ± 0.004 ± 0.005 

= 0.020 ± 0.008 ± 0.006 

0.066 ± 0.009 ± 0.008 
= 0.042 ± 0.007 ± 0.006 


5.0 ± 0.7 ± 0.8 

0.49 ± 0.28 ± 0.36 
(< 1.2 at 90% C.L.) 
0.39 ± 0.16 ± 0.12 
(< 0.7 at 90% C.L.) 
5.6 ± 0.8 ± 1.5 
3.5 ± 0.6 ± 0.9 


Si B(B»^K+tt-) 
f s B(.B°^n+n-) 


Si B(B»^K+,r-) 
B(B°^K+ K~) 


B(S U ^K+7T-) 
Sa B(A°^ P K-) 


Si B(B«^K+ir-) 
Sa S(A^pTr-) 


Si B{B"^K+tt~) 



consistent with the previous upper limit (< 5.6 x 10~ 6 at 90% C.L.) based on a subsample 
of the current data [5], and agrees with the prediction in Rcf. [15], but it is lower than 
most other predictions [12, 45, 46]. 

The -B° — > 7r + 7r~ upper limit improves and supersedes the previous best limit [5]. 
The present measurement of B(B° — > K + K~) is in agreement with other existing mea- 
surements and has a similar resolution [43], but the resulting upper limit is weaker due 
to the observed central value. The sensitivity to both B° — > K + K~ and B® — > n + TT~ is 
now close to the upper end of the theoretically expected range [12, 13, 14, 45, 47]. 

We also report the first branching fraction measurements of charmless decays. 
They are significantly lower than the previous upper limit of 2.3 x 10~ 5 [50], and in rea- 
sonable agreement with predictions [17], thus excluding the possibility of large (O(10 2 )) 
enhancements from R-parity violating supersymmetric scenarios [51]. Their ratio can 
be determined directly from our data with greater accuracy than the individual values. 
For this purpose, the additional pr > 6 GeV/c requirement is not necessary, and we can 
exploit the full sample size, obtaining B(A° b -> P TT-)/B(A° b -> pK~) = 0.66 ± 0.14 ± 0.08, 
in good agreement with the predicted range 0.60-0.62 [17], but in disagreement with the 
recent prediction in Ref. [52]. 

The dominant systematic uncertainties of all measurements presented here are due 
to finite size of control samples and are expected to reduce with future extensions of the 
measurements. 

In summary, we have searched for rare charmless decay modes of neutral 6-hadrons 
into pairs of charged hadrons in CDF data. We report the first observation of the modes 

— > K~tt + , A° — > p7r~, and A° — > pK~ , and measure their relative branching fractions. 
We set upper limits on the unobserved modes B° — > K + K~ and B° s -4- 7r + 7r~ . 
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